DLM在WinBUGS中的实现

贝叶斯动态线性模型中的  linear growth model，本要编了一个小程序，想实现一下，结果出错了，如下，那位高手帮我解答一下吧！
model{
   for (t in 2:n){
      y[t]~dnorm(mean.y[t],v)
mean.y[t]<-mean.mu[t]
      mu[t]~dnorm(mean.mu[t],w1)
mean.mu[t]<-mean.mu[t-1]+mean.beta[t-1]
      beta[t]~dnorm(mean.beta[t],w2)
      mean.beta[t]<-mean.beta[t-1]
}
    v~dnorm(mu0,tau0)
     w1~dnorm(mu0,tau0)
     w2~dnorm(mu0,tau0)
       mu0~dnorm(0,1.0E-6)
      tau0~dgamma(0.001,0.0001)
for (j in 2:m){
      pred.y[j]~dnorm(pred.mean.y[j],v)
pred.mean.y[j]<-pred.mean.mu[j]
      pred.mu[j]~dnorm(pred.mean.mu[j],w1)
pred.mean.mu[j]<-pred.mean.mu[j-1]+pred.mean.beta[j-1]
      pred.beta[j]~dnorm(pred.mean.beta[j-1],w2)
      pred.mean.beta[j]<-pred.mean.beta[j-1]}
}
#data
list(n=35,y=c(0.00,-0.70,-1.85,-3.00,-4.60,-5.55,-6.55,-13.90,-16.90,-19.25,
-21.60,-23.30,-24.85,-28.80,-33.35,-35.40,-37.10,-36.90,-36.80,-37.40,-37.70,-37.60,-37.60,-38.40,-38.50,
-38.00,-38.45,-39.05,-39.65,-40.70,-40.45,-40.65,-40.85,-40.70,-40.75),m=40,pred.mean.beta=c(1))
#inite
list(mu=1,mean.beta=1,v=0.1,w1=0.1,w2=0.1,pred.y=c(NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA
,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA)

咱没有人理我啊，快急死了

model{
   for (i in 1:n){
      y[i]~dnorm(mean.y[i],v)
    mu[i]~dnorm(mean.mu[i],w1)
     beta[i]~dnorm(mean.beta[i],w2)
mean.y[i]<-mean.mu[i-1]+mean.beta[i-1]
      mean.beta[i]<-mean.beta[i-1]
         }
       v~dnorm(mu1,tau1)
     w1~dnorm(mu2,tau2)
     w2~dnorm(mu3,tau3)
for ( k in 1:3){
      mu[k]~dnorm(0,1.0E-6)
      tau[k]~dgamma(0.001,0.0001)}
}
#data
list(n=35,y=c(0.00,-0.70,-1.85,-3.00,-4.60,-5.55,-6.55,-13.90,-16.90,-19.25,-21.60,-23.30,-24.85,-28.80,-33.35,-35.40,-37.10,-36.90,
-36.80,-37.40,-37.70,-37.60,-37.60,-38.40,-38.50,-38.00,-38.45,-39.05,-39.65,-40.70,-40.45,-40.65,-40.85,-40.70,-40.75))

#inite
list(mu=1,mean.beta=c(1),mean.mu=c(1),v=0.1,w1=0.1,w2=0.1)

出现错误mean.mu没有定义或在数据中没有设置

看来没有人理我了